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A differential ranging location system is de- 
scribed which uses a modified timc-of-anival tech- 
nique to detennine the location of a frequency- 
hopped spread spectrum radio signal. The trans- 
mitter simultaneously transmits two radio frequency 
carriers having different frequencies such that a 
phase difference is observed between the two carri- 
ers at a distance from the transmitter. The phase dif- 
ference is proportional to the range from the trans- 
mitter that the carrier signals have observed. The 
tw carrier signals from the single transmitter are 
received by at least three and in special cases four 
base stations which calculate the differential time 
of arrival based on the phase differences of the re- 
ceived carrien. The calculated phase differences 
are then sent to a central location which locates the 
position of the transmitter tmsed upon a planer hy- 
pertx>lic location algorithm. Since a large number 
of narr w band frequencies across a wide spectrum 
are used in a frequency hopping transmission pro- 
tocol, the accuracy in determining the phase differ- 
ences between the two carriers is increased and the 
immunity to interference is also increased. 
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DIFFERENTIAL RANGING FOR A FREQUENCY-HOPPED 
REMOTE POSITION DETERMINATION SYSTEM 

FIELD OF THF INVFMTIOki 

The present paiem application relates to radio communication 
and transmitter location system. 

BACKGROUND OF THF. INVFNTION 
Location of radio transmitters is a technology that has been 
used for many years and in many applications. For radio location, the use of 
mobile transmitters and a plurality of fixed base station receivers is generally 
known in the art where the fixed receivers locate the source of the radio signal 
by triangulation of the radio signals. Triangulation is a technique where a 
plurality of fixed base station receivers each determine the approximate 
15 direction of the radio signal from the transminer. The approximate direction 
lines are then drawn on a map and the location of the transmitter is defined at 
the cross point of these lines. Since the direction lines arc approximate, and 
typically three fixed base station receivers are used, the intersecting area 
between the three lines is usually shaped like a triangle since the lines rarely 
20 meet at a single point. The transmitter is then located within the triangle. 

Triangulation of radio signals is a costly undertaking since at 
least two receiving antenna stations must be situated at different diverse 
geographic locations to perform even rudimentan.' triangulation. Effective 
triangulation is more realistically achieved with three receiving locations and 
25 accuracy is increased even more with four receiving towers. The more 
receiving stations constructed, the better the accuracy in locating the 
transmitter, but also the great the cost. 

One such direction finding technique is described in US Patent 
No. 5.379,047 entitled "REMOTE POSITION DETERMINATION SYSTEM" 
30 issued January 3, 1995. In that patent, the direction of the incoming radio 
signal is determined by using an array of receiving antennas to measure the 
phase difference between the antennas of the array and thereby calculate the 
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incident dircclion of the radio carrier signal. Althouuh lliis patent has many 
novel and uscliil advantages, it requires that each base station use an array of 
antennas and at least two and typically three bases. Quite olten the cost ol' 
building and maintaining these base stations can be quite high, especially if 
they are located in urban areas where roof-top space must be leased and the 
cost of the lease is determined by the number of antennas placed thereon. 

Another disadvantage of iriangulation lor location of radio 
transmitters is the construction of such a system in sparsely-populated areas. 
In the example of a highway across a desert area, it is cost prohibitive to 
install a large number of base stations to track a narrow stretch of road. The 
base stations would have to be placed on both sides of the road and spaced 
quite far back from the road. 

There is a need in the art to locate a radio transmitter using at 
least three and in special cases four one-antenna receivers. A receiver with 
15 one antenna will minimize the construction and lease costs of a location 
system. There is a further need in the art to determine the location of low 
power, frequency hopped spread spectrum radio transmitter. 

SUMMARY OF THF lMVFMTir>Ki 
-0 The present invention solves the above-mentioned problems 

with locating transmitters and other problems which will be understood by 
those skilled in the art upon reading and understanding the present 
specificalion. The present invention uses a modified time-of-arrival technique 
to determine the location of a frequency hopped spread spectrum radio signal, 

25 The transmitter simultaneously transmits two radio frequency carriers having 
different frequencies such that a phase difference is observed between the two 
carriers at a distance from the transmitter. The phase difference is 
proportional to the range from the transmitter that the carrier signals are 
observed. The two carrier signals from the single transmitter are received by 

30 at least three and in special cases four base stations which calculate the 
differential time of arrival based on the phase differences of the received 
carriers. The calculated phase differences are then sent to a central location 
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which locales the posilion of liic transmitter based upon a planer hyperbolic 
location algorithm. Since a large number of narrow band frequencies across a 
wide spectrum arc used in a frequency hopping transmission protocol, tlic 
accuracy in determining the phase differences between the two carriers is 
increased and the immunity to intcrlerence is also increased. 

BRIEF DESCRIPTION OP THF r>n a 

In the drawings, where like numerals refer to like components 
throughout the several views. 

Figure 1 is a diagram of an integrated frequency hopping spread 
spectrum communication and location system; 

Figure 2 is a block diagram of a reverse paging and mobile 
location system using the infrastructure of a standard paging system; 

Figure 3 describes the synchronization and message format of 
1 5 the outgoing paging signals from the base stations; 

Figure 4 describes the formal of the frequency-hopped spread 
spectrum signal transmitted by the remote mobile units; 

Figure 5 is a hyperbolic plot of the differemial ranging of the 
present location system; 

2^ 6 is a diagram of the differential time-of-arrival aspects 

of the present location system; 

Figure 7 shows the computations used to perform location using 
differential time-of-arrival ranging; 

Figure 8 describes in more detail the principle of transmitting 
25 two simultaneous frequencies to two receiving base stations- 
Figure 9 is block diagram of the differential ranging transmitter: 
Figure 10 is a diagram of the IF circuit of the differential 
ranging receiver in the base station which uses two separate converters: and 

Figure 11 is a phase plot showing the effect of sources of errors 
30 on the differential ranging technique. 
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i)i:taiij:i) di-scrii'Tion of rm: im^.i.-h krh,) hmhoi^imwk.t 

In tiK- followiiiL. cleciilcci descripiion of ihe preferred 
cml-KHlimenl. reference is made u> ihc acccMnpanyinij drawings wiiieli lorm n 
pan hereof, and in which is shown hy way of illustration specilk preferred 
5 embodiments in which the inventions may be practiced. These cmbodimenls 
are described in sufllcienl detail to enable those skilled in the an to practice 
liic invcmion. and it is to be understood that other embodiments may be 
utilized and thai siruciural. logical and electrical chances may be made 
without departing from the scope of the presem inventions. The 

10 following detailed description is, therefore, not to be taken in a limiting sense, 
and the scope of the present inventions is defined by the appended 
claims. 

System Over\^iew 
In the preferred embodiment of the present invention, the base 
1 5 stations and central site g^e similar to the base stations 

and central site described in U.S. Patent Number 5,430,759 entitled "LOW- 
POWER FREQUENCY-HOPPED SPREAD SPECTRUM REVERSE 
PAGING SYSTEM" issued July 4, 1 995. and the processing .structure of this 
sy.stem operates similarly to the present system. The prci^cnl invention allows 
20 the accurate location of the frequency hopped, spread spectrum transmitters 
even in noisy and multi-path environments such as those found in urban 
environments and irregular terrain environments. 

The present invention is directed to a location system for use in 
locating remote mobile units (RMU's). The remote mobile units may be 
25 located in motor vehicles, located on the person of people, carried within 
containers or packages, or any number or variety of mobile carriers. The 
location system comprises a plurality of base stations (at least three and in 
special cases four) which can transmit signals to the remote mobile units and 
receive signals from the remote mobile units. The base stations may send 
30 messages and information to the remote mobile units, the base stations may 
send alert or interrogation commands to activate the remote mobile units, or 
the remote mobile units may be activated locally without any prompting. 
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Thus the remoie mobile units may he activated locally (for example by the 
wearer or vehicle) or remotely by the base station. The remote mobile units 
may be used by a person as a reverse pager, as an emergency locator or as a 
communication device. The remote mobile units may also be used as part of 
a vehicle to locale the vehicle ii' stolen, lor vehicle tracking or as an integral 
communication device. The remoie mobile unit may also be activated in the 
case of an accident by a deceleration switch. ABS system, air-bag or other 
crash-activated switch to automatically summon help. 

The remote mobile units receive messages as a standard paging 

10 device over licensed airwaves using a standard paging infrastructure. The 
remote mobile units, when activated, transmit low-pov^-er (less than one watt), 
frequency-hopped, spread -spectrum communication signals. The transmitted 
signals from the remote mobile units are received by the base stations which 
are then used to locate the remote mobile units. The base stations are 

15 equipped with a single receiving antenna. The location of the remote mobile 
units is determined by differential ranging using phase difference comparisons 
between two carrier signals simultaneously transmitted by the remote mobile 
unit. Unique algorithms are performed on the received signals at the base 
stations to eliminate the ambiguity caused by signal interference, phase 

20 ambiguity caused by errors and noise and multi-path rettections. 

The frequency hopping transmitter of the remote mobile unit 
produces two sine waves for each frequency hop. The phase difference 
between the two sine waves will vary with distance from the transmitter which 
is proportional to the range from the transmitter. There is an ambiguity factor 

25 in the range determination due to the phase differences between the two sine 
waves. The ambiguity is minimized by transmitting at least 53 hops per 
transmitted message so that there is a different phase difference within each 
frequency hop. Each hop change produces a step in frequency such that 53 
steps of various step distances in frequency are produced. 

30 At the receiver, the frequency hopping carriers are 

downconverted. sampled and digitized for processing in the digital domain. A 
range calculation is determined by comparing the phase differences of the two 
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sine waves of ihc two frequencies within each hop. The ambiguiiics are 
minimized by comparing the phase difTcrcnces across the 53 hops in 
synchrony with the iiopping protocol. Tlie difTerential range between two 
bases is determined by collecting the phase dilTercnccs Irom all the hops in 
two bases and pcrlbrming in the central station an Inverse Fourier Transform 
(IFT) which will produce a single value Ibr the differential range beiwee.i the 
two bases to the transinilier. 

Each base station is equipped with a very accurate rubidium 
atomic clock so that the exact time of arrival of the frequency hopping signals 
is determined with reference to a very accurate standard time. Those skilled 
in the art will readily recognize that any highly-accurate clock means may be 
substituted for the preferred rubidium clock such as a GPS receiver. The 
central station compares the time of arrivals with t..e knowledge that each 
base station has made its calculations based on the same time standard. The 
1 5 lime of arrival information is then used to plot possible locations of the 
transmitter using a hyperbolic plot. The comparison of the lime of arrivals 
from two base stations will produce a hyperbolic plot where the base stations 
are the foci. By plottmg the time of arrival numbers on the hyperbolic plot 
and crossing the hyperbolic plots between two or more base stations, the exact 
20 location of the iransniilter can be determined. 

To assist in calibrating the differential 'ranging system with the 
local base stations, a transmitter having a known location is used to 
periodically transmit a signal which is then used to synchronize the base 
stations. Tlie calibration signals is used to adjust the local atomic clocks of 
25 the base stations. It is possible to use less accurate clocks in the present 
system by increasing the number of calibrations performed in a day. In a 
calibrated system the accuracy of location has been measured at 26 meters at 
minimum received signal strength of -138dbm when operating within the 
parameters described below. 

'^^^ problems associated with multipath rellections and fading 
are known in the art. The multi-path problems and fading problems are 
usually frequency dependem so the use of a large number offrequency hops 
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N irlunlly climinalcs .hcsc problems. The ambigimy problems arc also soK-cci 
in n lashion .similar lo Jhat ileseribecl in the aroremenlioned coponclin,.. palcnl 
application. 

Summary of Hyperbolic Location with nin ercmial Tim..- of Arrival 

I hc operation of the present differenlial ranging system may 
also be termed differential time-of-arrival (TOA) since the differences in ihc 
arrival times at the base stations are used to determine the position of the 
receiver. As is described in more detail below, the TOA of the lv,o carrier 
frequency signals for each remote mobile unit are estimated at each base 
station for each hop. The TOA is proportional to the range R of the signals. 

The operation of the present system is basically performed in 
two major steps. In the first step, estimated range is determined for each 
base station BS„, where n is the base station number. Due to the inaccuracy 
15 of the clock in the transmitter, the range by itself is not an accurate indicator 
of the distance of the remote mobile unit from the transmitter. Since the base 
stations are equipped with very accurate clocks which are periodically 
synchronized to a common time base, the base stations can report the exact 
time of arrival of the signals. 

In the second major .step, the time of arrival information from 
each base station and the range information (calculated across the multiple 
hops using an Inverse Fourier Transform) are sent to a central station where a 
range differential between two base stations is plotted, as shown in Figure 5. 
The range differemial forms a hyperbolic plot, for example 501a and 501b, 
between base station BS, and base station BS,. We know that the range R, is 
greater than the range R, so that the hyperbolic plot 501b in dashed lines in 
Figure 5 can be ignored. A second hyperbolic plot 502a and 502b is formed 
between base station BS, and base station BS,. The intersection point 
between the two hyperbolic plots is the location of the remote mobile unit 
30 RMU. 



20 



25 
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Comparison of Dirreronlial Uanuinu u> (il>S 
The operation ol^ ihc present dilTcrcntial raniiinii system is quite 
ililVereni IVom the C'iPS (C}lobal Posiiioninu SaieUile) system. I lie CiPS is a 
naviiiaiion system that provides location for the vehicle that carries it. In 
5 order to monitor the vehicles, the system requires a communication system for 
delivering the location to a central ofllce. The navii-ation process of the GPS 
is virtually the inverse process of the dual frequency dilTerenlial ranging 
system of the present invention. The theoretical accuracy of the present 
invention is only dependent on the bandwidth, signal to noise ratio and some 
10 ambiguities described below. Although commercial GPS systems are made 
inaccurate on purpose, a comparison can still be made. 

There are several facts in favor of the present differential 
ranging system for the task of remote mobile unit location over GPS. The 
GPS system requires direct line of sight to the satellites, since it has a very 
1 5 sensitive receiver. Differential ranging is transmitting in the 900 MHz band 
and thus does not require line of sight: it could operate from within buildings. 
Differential ranging does not require any other device other than the low cost 
RMU on the vehicle which makes it significantly lower cost. The power 
consumption of the differential ranging system is signincanily lower since it 

20 requires only a low power transceiver, thus it is. operates on batteries for a 
long time (over 1 month on commercial batteries). The differential ranging 
system operates using frequency hopping and is less susceptible to multipatli 
and interference than the GPS which operates in direct sequence spread 
spectrum. From measurements performed on the present direction finding 

25 system, it is by far less susceptible to ground multipath than the GPS. The 
complex receivers of the system are only in the base stations. The 

complex receivers of the GPS are in every vehicle. 

Frequency Hoppinp System Ovprvi^.^, 
30 In the preferred embodiment of the present invention, 

transmitters and base stations similar to those used in the present invention are 
described in the aforementioned copending U.S. Patent Number 5,430 759 
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cmiiicci "Low-powi-K i-Ki:oiJi:NCY-ii()i>i»i;n si>ui-:ai) si>i:(trum 

in-VI-KSI- PAGINfJ SYS Tl-M" issued .h.ly 4. The infrastmclurc of 

lliis reverse pjitiint' syslcm serves :vs the hnsis ol" ilie present :»nd the 

reverse pagers of this syslcm operate similarly in tlie present svstcm using a 
5 single carrier frequency instead of ihe two carrier signals of the present 

invention. The system «l'«ws the accurate location of the reverse 

pagers even in noisy and multi-path environments such as those found in 
urban environments and irregular terrain environments. 

The environmenl in which the present system is operated is 
10 shown in Figure 1. The system is designed to operate in a 

communications environment known as NexNet™ in which a plurality of base 
stations service a particular area using the NexNet™ communications system. 
Each remote mobile unit receives standard paging messages as a downlink and 
can provide talk back messages as an uplink channel using frequency hopping 
1 5 spread spectrum radio frequency communication. The paging transmitter of 
the base stations 113 typically transmits with 200 Watts of power to send 
messages using a standard paging format such as POCSAG or other industry 
standard formats. Each remote mobile unit contains a paging receiver with a 
sensitivity of -115 dBm. a frequency-hopping dual frequency transmitter with 
20 a transmit power of one Watt and double side band modulator: Each base 
station includes a receiving antenna for the frequency hopped carrier signals 
and special processors for determining the phase difference between the two 
carriers of the double side band modulation. 

Figure 2 depicts the major components of the two-way paging 
25 system in the aforementioned U.S. Patent Number 5,430,759 entitled "LOW- 
POWER FREQUENCY-HOPPED SPREAD SPECTRUM REVERSE 
PAGING SYSTEM" issued July 4, 1995. The reverse paging terminal 1 10 at 
the central site operates to provide synchronization and messaging information 
through the paging terminal 1 16 to the reverse pagers 100 (also known as 
remote mobile units 100) via direct links to the base stations BS„ BSj and 
BS, through ground based radio links (not shown) or through a satellite 
uplink/downlink using a geostationary satellite (not shown). The base stations 
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I^S,. MS, and BS, include iransmil and receive lowers 11 3a. 1 13b and 1 1 3c. 
respectively and base station terminals 201a. 201b and 201 c. respectively. 
Terminals 20la-201c arc required for producing ihe accurate .synchronization . 
inlbrmalion needed to be transmitted to the remote mobile units and for local 
processing of the received messages for transmitter location. This 
synchronization information is used to coordinate the iVequency hopping 
transmissions and to coordinate the response of messaging from the plurality 
remote mobile units 100 so as to minimize collisions within groups of remote 
mobile units and eliminate collisions between groups of remote mobile units. 

10 

Base to Remote Synrhr onization 
Standard paging messages sent from the base stations BS,. BS, 
and BS, to the plurality of remote mobile units 100 are. in the preferred 
embodiment, sent as digital data encoded in the POCSAG paging standard. 
15 These messages may be used to interrogate the remote mobile units to activate 
the, remote mobile unit lo allow the base stations to begin the location process. 
Typically the paging channel has a center. frequency of 143.160 MHz. with an 
NRZ FSK data rate of 5 1 2 bps or 1 200 bps. Other bit rates such as 2400 
baud (bps) are also feasible. Figure 3 describes the POCSAG paging 
20 communications protocol as modified for use by the preferred embodiments of 
the present invention. In the top line of Figure 3. a greatly compressed time 
line of digital data transmitted according to the POCSAG protocol is shown. 
Batches of messages are transmitted in groups as shown in the details in the 
subsequent lines below the top line of Figure 3. In the second line of Figure 
25 3. a 1.0625 second interval (for 512 baud) is shown in which 544 bits are 
transmitted as a single batch. The batch is preceded by a synchronization 
code word SC as shown in the third line of Fiaure 3. 

The synchronization code word within each batch is followed 
by eight frames of digital data. Each frame is divided into two portions, an 
30 address portion and a message portion. The address code word of the 

message of frame 2 of Figure 3 is shown in line 4 while the message code 
word of the second half of frame 2 is shown in line 5. The address code 
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word is preceded by a digital zero followed by i X address biis. two lunction 
bits and 10 check bits. 1 be address code word is Collovvcd by an even parity 
bit. The message code word portion of the frame is preceded b> a digital one 
followed by 20 message bits which are followed by 10 check bits and a single 
5 even parity bit. fhus each frame is comprised of 64 bits divided into two 32 
bit sections. 

Synchronization of the base station terminal 200 and the remote 
mobile units 100 is necessary to ensure the units 100 are transmitting at the 
same time that the base stations are listening. Synchronization is also 
10 necessary to coordinate the division of the large number of remote mobile 
units into groups so that members of one group use different frequency 
hopping patterns from members of other groups. Synchronization of the 
remote mobile units 100 is accomplished by inserting a special frame into the 
POCSAG data which is used to synchronize the units. 

The purpose of synchronization between the reverse paging 
terminal 1 16 and the remote mobile units 100 is to determine where along the 
pseudo random noise code the frequency hops are to be followed and lo 
determine the exact times for transmitting frequencies from within any of the 
hops. This also enables the dynamic changing of a remote mobile units group 
membership such that if one group is experiencing a large number of 
collisions due to simultaneous transmissions, the reverse paging terminal 110 
may re-allocate some of the remote mobile units within that group to new 
groups to minimize collisions. 

Referring once again to Figure 3. eight frames of information 
25 are transmitted in each burst using the POCSAG formal. Remote mobile units 
100 may be assigned to a specific frame within the transmission so that the 
remote mobile units, once recognizing the synchronization code word, can 
scan a specific frame for that remote mobile unifs address. Once the address 
is found, the remote mobile unit can determine any group changes that may be 
30 required to re-allocate that remoie mobile unit to a different group. In 

addition, the POCSAG format is used to transmit a fine time synchronization 
code. The fine synchronization code is a transmission of a time pulse at an 
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exact lime synchronized lo a CJ!>S (Cilobal Positioning System) clock to 
synchronize all the remote mobile units 100 for time of transmission, l or 
example, periodically during the day the reverse paging terminal xvill send a 
synchronization code within the POCSAG code word which is sent 
precise time. In order to ensure that a precise time pulse is sent, the rex. 
paging terminal 1 10 receives accurate time information using a GPS antenna 
to receive accurate time of day information. The time used to send the 
synchronization pulse is when the day clock reaches exactly some multiple of 

0.9 seconds in the preferred embodiment in tu:^ „ u • • . 

H nooaimeni. in this synchronization information. 

20 bits of information are transmitted to give the accurate time of da>- 

information. 

In each of the remote mobile units 100, the microprocessor 
compares this accurate time pulse which will indicate the exact time of dav 
and compare it to its own day clock. The clock within each microprocessor is 
1 5 accurate down to a few milliseconds, but the time at which the 

synchronization pulse occurs should have a resolution much fmer than that 
such as down to 0.1 milliseconds for time of day. In this fashion, each of the 
microprocessors in each of the remote reverse paging devices can periodically 
realign its day clock to know within a millisecond the exact time. Each ' 
microprocessor does not actually realign its clock but changes a clock offset 
within memory so that it understands how far off its own internal clock is and 
can make the adjustment when using that clock to determine when to start 
transmitting information. 

The synchronization pulse is only transmitted everv tew 
minutes. However, the resolution of the start of the message indicatinu the 
synchronization pulse is very accurate, it being transmitted at 0.090000 
seconds GPS time after a fixed time of day. such as 12:00 GMT. This GPS 
time is accurate lo at least within 1 00 nanoseconds. 

An overview of the transmission format of the remote mobile 
unit is shown in Figure 4. The actual transmission of information from the 
remote mobile units 100 is done using Differential Bi-Phase Shift Keving 
(DBPSK) modulation on a frequency hopped carrier of less than one wau. 
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The transmission of inlormaiion Ironi liic rcmoic mobile units KK) on liic 
frequency hopped carrier may also be done using I-requcncy Shift Keyinu 
(1-SK) modulation. 'I he same inlormaiion is modulated on the two carriers in 
the preferred embodiment of the present invention allowing the base station to 
5 decode the information from either carrier or both carriers. The dual carriers 
arc preferably two different narrow hop frequencies. 

Typically a single transmission consists of 53 hops or 53 
changed frequencies selected from a list of narrow band frequencies. The 
frequency selection is based on a pseudo-random noise code list pointing to 

10 the frequency selection list. The synchronization information tells the remote 
mobile unit 100 where along the pseudo random noise code it should be 
synchronized for transmission of its message and tells exactly the time of day 
so that the remote mobile unit 100 knows exactly when to start transmitting 
the specific frequency so that the base stations BSl. BS2 and BS3. labeled 

15. 1 1 3a, 1 1 3b and 1 1 3c. respectively, are looking for that frequency at the same 
lime. 

In operation. 200 frequencies -are used by the remote mobile 
unit 100 and the base stations and internajly. stored in a list numbered Fl 
through F200. For a specific message. 53 frequencies will be used to transmit 

20 the entire message. These 53 frequencies are selected based on a 1.000 
member pseudo-random noise code. 

The use of the accurate synchronization signal periodically 
broadcast via the outbound paging signal enables the remote mobile units to 
use lower accuracy components thus reducing the manufacturing cost of 

25 remote mobile units. For example, high accuracy crystals to track the time of 
day within the microprocessor are available with an accuracy of three pans 
per million. Thus, a time drift of approximately three micro seconds per 
second or 180 microseconds in a minute is the known drift. There are also 
time inaccuracies which are introduced due to variable path length between 

30 the paging tower to the remote mobile unit. By employing crystals which are 
cheaper and have an accuracy of the order 50 parts per million, the amoum of 
time-of-day drift normally wouldn't be tolerable. However, by using the 
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synchronizalion information iransmiiicd on a regular basis from ihc reverse 
paging terminal, the microprocessor can continually correct its own internal 
day clock so that accurate time of day measurements are always mamiained. 
,The microprocessor estimates the momentary inaccuracy of the cry stal b\ 

5 tracking the drift across several synchronizalion transmissions and dynamicalK 
adjusts for the frequency drift of the crystal and the offset using internal olTset 
registers for accurate time of day information. 

Counters are employed within each microprocessor of the 
reverse paging units to compensate for the offset of the frequency based on 

10 the synchronization time information. There are generally two major factors 
which affect the drift in a crystal: temperature and acceleration. Most of the 
drift is due to temperature, and the remaining drift components are negligible. 
The frequency drift in a cr>'stal due to temperature is very slow, on the order 
of 50 Hz over 10 seconds. During a single, day the temperature can change by 

15 20 or 30 degrees fahrenheit. requiring a lime update from the GPS clock 
approximately every five minutes. 

Remote Mobile Unit Transmisginn Formal 
The signal sent from the remote mobile unii 100 to ibe..base 
20 stations is a spread-spectrum, frequency-hopped transmission using differential 
bi-phase shift keying (DBPSK) modulation on the frequency -hopped carrier to 
transmit digital information. The transmission of information from the remote 
mobile units 100 on the frequency hopped carrier may also be done using 
Frequency Shift Keying (FSK) modulation. The frequency hops are relatively 
25 slow, the frequencies transmitted are very narrow and the transmission power 
is extremely small. 

The maximum peak output power of transmission from remote 
mobile unit 100 is limited to less than one Watt to allow use of the 902-928 
MHz ISM band in the United States without the need for licensing the remote 
30 mobile units as allowed by FCC regulations defmed in 47 C.F.R. §15.247. 
Those skilled in the art will readily recognize that other frequency bands and 



wo 96/25673 



PCT/GB96/00270 



1 



10 



iransmissions power levels may he cmployetl iJepeiulini. upon IVC licensing 
iciniiroiiients or olhcr frequency licensinii rcqiiircnienis of other nniionaliiies. 

I he use of mi aeeuraio crystal u. eontrol each iVcMuency of 
iransmission is required within each remote niohiie unit 100. l or example, 
higli accuracy crystals to transmit the narrow bandwidth frequencies used for 
the frequency hopped transmissions are available with an accuracy of three 
pans per million. At 900 MHz. a 3 ppm drift would place a single frequency 
somewhere within a 2.7 KHz band. To tolerate frequency drift due to aging 
and temperature, each individual frequency of the frequency hopped signal is 
allocated to a 7.5 KHz band or channel, even though the actual frequency is 
on the order of 200 Hz wide skirt within this 7.5 KHz allocated bandwidth. 
Those skilled in the art will readily recognize that by using alternate 
components, the frequency channels (individual frequency of the frequency 
hopped signals) of 7.5 KHz allocated bandwidth may be wider or more 
15 narrow depending upon the overall allocated bandwidth for the system. For 
example, 1 KHz or less bands may alternatively be allocated per channel. 

Tests on this system have shown that by processing the 
received signals at the base stations entirely in the digital domain using the 
conibination of unique Fast Fourier Transform algorithms 

locate and retrieve the frequency hops and by using a 
combination of unique confidence algorithms with a plurality of error 
correction codes, the receiving base station is able to pull the response 
information from a very low power signal from a distance of up to 45 
kilometers (28 miles) in a flat terrain. In a rather noisy urban environment, a 
25 range of 24 kilometers (15 miles) is the norm. The information within the 
signals is accurately decoded even in severe multipath and noise conditions. 

As shown in Table 1. the remote mobile unit message format 
consists of a preamble and the message body spanning a total of 53 frequency 
hops. Those skilled in the art will readily recognize that longer messages may 
30 be transmitted using the preferred embodiment of the present invention, and 
the messages format described here is illustrative but not limiting. Much 
longer message hops to transfer more digital data is also implemented but not 
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clcsciihcci here. Of course, those skilled in iIk- :irl rciidilv reeoiriii/.e iIkiI 
shorter n-iess:mes than those dcscrihcd helow :iie e<iiiall\ pcissililc lor tlie 
|irclerred embodiments of the present invention. I lu- niossaiio Icnetli nnd 
ninnber ol" transmission hops arc a matter of desitiii choiee. 
5 l lie message preamble consists of a predelined eode of ones 

and zeros to gel the attention of the base unit receiver to pull the mcssauc out 
of the noise. The preamble consists of 165 bits transmitted across 5 hops, that 
is. transmitted using DBPSK (Differential Bi-Phasc Shift Keying) or 
Frequency Shift Keying (FSK) on five different frequencies selected from the 
10 frequency list with the specific frequencies selected based on the PN (Pseudo- 
random Noise) Code list stored within the remote mobile unit. The sequence 
location within the PN code that the remote mobile unit will begin to follow 
is based on the synchronized time of day. Within a single hop (a single 
carrier frequency), the carrier phase is modulated 33 times to encode the 
15 predefined one-zero pattern of the preamble. 

The message body follows the preamble and consists of three 
groups of data. Each group consists of 30 actual data bits so that the entire 
message is, in the preferred embodiment of the present invention. 90 total data 
bits (although other bit length: messages may be chosen). The actual data 
20 encoded within these 90 bits is described above and may be in any convenient 
coded format. Those skilled in the art will readily recognize that a wide 
variety of message formats and encoding of the data bits may be used without 
departing from the scope of the present invention. The encoding 

described here, however, has been proven effective in retrieving the data bits 
25 buried in background noise with a high degree of accuracy and a low actual 
error rate. 

TABLE 1: Remotie Mobile Unit Message Format: 
Preatnble is 165 bits (33 bits x 5 hops) 
Message is 48*33 transmitted bits 
30 (Message is 90 bits actual data) 

[] = One Frequency Hop 

I ' > ' , ' I ' ' ' ''II — I I I — 1 — I — \ — I — I — I — \ r- 

preamble I Message Body Spread Over 48 Hops 
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l:ach of the three groups ol' message daui (?>() bits eaeh) arc 
new encoded using a standard 30.63 BCII eode and with a single parity bii 
added to Ibrm a 64-bil word. This encoding decreases the error rale iVom 10 - 
5 lo 10 \ This encoding, documented and understood by those skilled in the an. 
can correct up lo 6 errors or delect up to 13 errors. Deieciion of corruption of 
a data word thai cannot be reconstructed will cause the base lo request a 
second transmission of ihe acknowleduemenl message. 

10 I nner Codinii and Interlea vino 

The inner coding of the message will protect the integrity of the 
message with an error rate as high as 25%. Each block of 64 bits of data 
(corresponding to a groups of 30 bits and earlier encoded by a standard 30.63 
BCH code) is split into two sub-blocks of 32 bits (sub-blocks A and B of 

1 5 Table 2). and a reference bit is added to each sub-block to assist the 
differential encoding to provide a reference bit to the DBPSK or FSK 
decoder. The 33 bit sub-blocks are transmitted over one frequency hop each 
and are replicated 8 times so that the 64.bit block traverses 16 frequency 
hops. In transmission, the 33 bit sub-blocks sire interleaved to further reduce 

20 loss of data, as shown in Table 3, where sub-blocks A and B of Table 2 

correspond to the first group of 30 bits, sub-blocks C and D. correspond to the 
second group of 30 bits, etc. The total message is 53 hops where each hop is 
180 msec in length making the duration of a single messace 9.54 seconds 
TABLE 2 : Interleaving Format for Sub-block 

25 

A = 1 reference bit and 32 data bits = 33 bits 
B = 1 reference bit and 32 data bits = 33 bits 



30 



[] = One Frequency Hop 

[Tbababababababab 
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TABLE 3: Inner Coding and Interleaving of Siab-blocks 

A = first 33 bits of isz block 
B second 33 bits of Isc block 
C = first 33 bits of 2nd block 
D 



E 
F 



= second 33 bit 



s of 2nd block 



= first 33 bits of 3rd block 
= second 33 bits of 3rd block 

Q = One Frequency Hop 



1 1 I 1 
preamble 
till 


A 


B 


A. . .B 


A 


B 


C 


D 


C . . . D 
X 


C 


D 


E 


f| . 


E 


F 


E 


F 



Those skilled in the an will readily recognize thai a wide 
variety of data interleaving may be utilized to effect better error tolerance and 
may be substituted for the interleaving described here. Such alternate 
substitute interleaving means are CIRC (Cross Interleaved Reed Solomon 
20 Code) used in CD (Compact Disc) recording media operating either at the 
block level or at the bit level. 



Sinule Hon Format 

The acknowledgment signals are transmitted by the remote 
25 mobile units 100 in a 1 .5 MHz band selected from within the 902-928 MHz 
spectrum. The 1.5 Mhz band is divided into 7.5 KHz channels to provide 200 
channels available in which the frequency hops can occur. Thus, each 
frequency hop is a channel 7.5 KHz wide in which a carrier frequency is 
transmitted. For example, channel one will have a frequency Fl at 
30 902.00375 MHz + -3.75 KHz, channel two will have its center carrier 
frequency at 902.01025 MHz + -3.75 KHz. etc. 

Each transmit frequency of each hop will thus be centered at 
the approximate mid-point of the assigned channel band: however, due to 
inaccuracies in the remote mobile unit circuits and reference crystals, the 
actual transmit frequencies will vary between units. If high quality crystals 
are used to accurately produce the required frequencies, very little drift off the 
center frequency will result. In the preferred embodiment of the present 
invention, low cost crystals are purposely employed to keep the per-unit 
manufacturing costs down. This will allow for a lowcr-ccst product sold to 



3> 
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ihe user which will increase market pcnclralion. Thus, reference cryslals are 
preferred which have a frequency accuracy of 3 ppm such thai al 000 MHz. 
the siaiislical drift would be approximately 2700 Hz. The cryslals center 
frequency within its nominal accuracy also drifts due to aging and temperature 
5 variations, but this drift is slow compared to the transmission limes so the 
drift during a single transmission due to these latter variants is unimportant. 

A single frequency hop is shown in Table 4. The 1 5 
millisecond guard time preceding each hop is primarily a sealing time for the 
oscillator circuits of the remote mobile units to allow the internal oscillator 
10 circuit to lock onto the new frequency between hops. Each hop is transmitted 
at a single frequency in which the phase of the carrier is either at 0 degrees 
phase or 1 80 degrees phase in reference to the phase of the reference bit 
immediately following the quiet or guard time. Thus the first bit is a phase 
reference bit followed by 32 data bits exhibiting either zero phase shift or 180 
15 degree phase shift to encode the data bits as DBPSK (Differential Bi-Phase 
Shift Keying). In an alternative implementation, each frequency hop may be 
modulated using Frequency Shift Keying (FSK) in which two frequencies are 
used to transmit data bits. One hop frequency may indicate a logical one 
while a second hop frequency may indicate a logical zero. The frequency 
20 shift is minor and the frequency differential is contained within a single hop 
channel. 

Each bit of DBPSK or FSK is a transmission of 5 milliseconds 
of the hop carrier frequency either in phase with the reference bit transmission 
or 1 80 degrees out of phase. 
25 TABLE 4 : Single Frequency Hop Format 

Guard time (quiet) = 15 ms 

Single Bit = 5 ms of carrier DBPSK / FSK 

33 Bits plus guard time = 180 ms 
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5ms 
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5ms 
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l-reciuencv 1 lonniiiL' SccuiencL> 
All of ihc remote mobile units in the market serviced hy ihc reverse 
paging terminal lor message or location finding use the same pseudo random noise 
code to determine the frequency hops, fhe pseudo random noise code is a digital 

5 code of 1000 unique numbers. In the preferred embodiment of the present 

invention, the pseudo random noise code is stored in memory of each of the remote 

. mobile units. Those skilled in the art will readily recognize, however, that a linear 
feedback shift register could be used to generate the pseudo random noise code on 
a real-time basis instead of using a look-up table which is presenth- in the preferred 

10 embodiment. 

The PN (pseudo-random noise) code list is stored in memory and 
maps to a frequency list. In the preferred embodiment of the present invention, the 
PN code list has 1,000 entries which repeat as a sequence. The control means of 
the reverse paging units continuously maintain a count of the proper location within 
15 this list. As described below, the time of day for all remote mobile units in the 
market served by the base terminal are periodically synchronized to ensure 

acknowledgment messages are synchronized to transmit the hop frequency at the 
proper time and to synchronize the location within the PN code list that each 
remote mobile unit will use to transmit. 

20 The 1,000 member PN code list maps to a 200 member frequency 

list. In order to allow a large number of remote mobile units to simultaneously 
operate in the same geographic market, the remote mobile units are divided imo 
groups and the groups are assigned different sequence segment locations in the 
-same 1.000 member PN list. Thus a remote mobile una from group one will begin 

25 transmitting a hop at a frequency determined from a first location with the PN 
code, while a remote mobile unit from group two may begin transmitting a hop at 
a frequency determined from a second location in the PN code. The remote mobile 
units from group one and group two will complete their respective 
acknowledgement messages in 53 hops. Preferably, the sequence of the PN code 

30 used to determine the frequencies of the 53 hops for the remote mobile unit of the 
first group will not overlap the sequence of the PN code used to determine the 
frequencies of the 53 hops for the remote mobile unit of the second group. More 
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preferable, the IVcqucncics chosen bused on the non-overlapping segments of the 
I^N code list are orthogonal such thai the same frequency is never used by two 
remote mobile units belonging lo different groups. 

In the preferred implementation, the LOGO member PN code list is 

5 divided into 160 hopping sequences. The remoic paging units are divided into 40 
groups with the members of each group synchronized to track the same location in 
the PN code list. The microcontroller of each remote mobile unit, regardless of its 
group membership, continuously runs through the repeating PN code sequence to 
stay in synchronization with the base unit and all other remote mobile units. Each 

10 group of remote mobile units is further divided into four subgroups such that the 
remote mobile units within each subgroup are assigned one sequence within the PN 
code list. Although the 53 hop sequence needed for each acknowledgement 
transmission may overlap the 53 hop sequence used by a remote mobile unit in 
another subgroup, the transmission sequences of a remote mobile unit of one group 

15 is chosen to not overlap the 53 hop sequence used by a remote mobile unit in 
another group. 

Base Station Desig n 

The analysis of the dual carriers and the decoding of the signals 
20 received by the base stations from the remote mobile units is done almost entirely 
in the digital domain. The carrier frequencies of the frequency hops are down- 
converted to a lower frequency in each base station and are then digitally sampled. 
The digital samples are then processed to locate the phase information of interest 
for direction finding and message decoding. Each base station is constructed with 
25 a plurality of digital signal processing based receivers which enable simultaneous 
message decoding and direction finding of a plurality of simultaneously 
transmitting remote mobile units. The dual carriers arc analyzed in parallel to 
determine the exact phase separation. Each bases station is also equipped with a 
rubidium atomic clock for accurate time-of-day calculation. Those skilled in the 
30 art will readily recognize that otlier, less accurate clocks such as GPS clocks can be 
used but the effective location resolution may be affected. 
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Dincrcnl i;il l ime ol' Arrival Svsit-m 
Figure 6 shows the basic principals of the Ditrcreniial Time of Arrival 
(DTOA) System. The example system is configured by at least four reccivin-: 
stations which operate at exactly the same clock. The transmitted siynal is received 
5 by the receiving stations al T,. 'I T,. and T, respectively. The base stations are 
located at (X Y (X Y (X Y ,) and (X Y J. The four time 
measurements from the four base stations are resolved in non ambiguous fasliion as 
follows: 

10 C = Velocity of Light. ^ ^ ^ 

X,. Y, are the coordinates of the transmitter, and are resolved directly from 
those 3 equations. T, is a dummy variable representing the time of transmission. 

(X,-X,)-^-h(Y ,-Y.)V=C-*(T,-T,)-\ (2) 
15 (X,-X.)= + ( Y,-Y.)- =C^MT,.T.V 

(X3-X.)^ + ( Y,- Y,)= =C^*(T,-T.)^ 
(X,-X.)^ + ( Y,- Y,)^ =C-\*(T,-T,)^ 

The four equations are dependent (rank 3). Therefore three linear equations could 
20 be from them (by subtracting row 1 from row 2. row 2 IVom row 3 and row 3 
from row 4): 

(beftnc: A P , , = P . - P, : where P = Oany parameter and CAT , , = A R . , ) (3 ) 
2AX X. + 2AY,3 Y. - 2CAR,, T. = x + Y,- - X - Y,- - C" (T,--T,^) (3a) 
25 2AX X . + 2AY,3 Y. - 2CAR,, T, = X + Y,- - X - Y,- - C- {T,--T,-) (3b) 
2AX,, X. + 2AY„ Y,-2CAR„T, = X,- + Y,--X,= - Y,- -C^(T^T/) (3c) 

Those three equations could be easily further reduced into two linear equations lor 
X , . Y , dependent only upon AR ij. 
-'0 By dividing equation (3a) by A R,, and equation (3b) by A R„ and 

subtract the results and do the same for equations (3b) and (3c) we get the 
following equations. 
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(Define: I).; = X ; + Y; - X - Y;) 

X, !A X,, /A R,, - A X,, /A R,,l + Y, [AY,, /A R,, - AY,, /AR,,! = 

{(D„- /AR„) - (D.J I AR,0 - AR,,]/2 
X, !A X,, /A R,. - A X,, /A R,,l + Y, lAY,, /A R,, - AY,, /A R,,! = 

{{D,,- /A R,, ) - ( D,/ / A R,,) . A R,,)/2 



The location could be resolved with only three bases, but it would 
result in an ambiguous solution which could be resolved by an a priori 
knowledge of the transmitter location. 

10 

Dual Frequency Differential Ranging Princip le 
The differential ranging system is based upon transmitting two 
simultaneous frequencies for each hop, with an alternating difference that solves the 
ambiguity. The principal is described in Figure 8. The mathematical derivation of 
15 the differential ranging is based upon the phase difference that is created when two 
frequencies with a small difference in their wavelength are propagating over 
distance between the transmitter and the receiver. The two signals are represented 



as: 



20 A, (t) = A„ ♦ Cos (wt + f) (5) . 

A, (ty= A„ * Cos ((w.Aw)t + <j)„) (6) 

The phase difference that is propagating through the distance R^ 
between the transmitter and the first receiver is: 

25 

Atp, = wR, /C + 4) - (w - Aw) R, /C - (|) 0 (7) 

The same applies also for the second receiver which is at a distance of R. 
from the transmitter: 

30 ' 

A ip: = wR. /C + <}> - (w - Aw)R, /C - ^ (8) 
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Solvini! the two equations: 

AK,: = R, - = (A(p, " A<p. )* C/Aw ^^)) 

5 A (P( and A tp. are measured by ihc digital receiver at each base siaiion. Aw 

is known a priori. The basic algorithm is as follows: 

1. Compute the phase difference of the two frequencies at the first base 
station. 

2. Compute the phase difference of the two frequencies at the second 
10 base station. 

3. Compute the differential range using formula (9) above. 

Ambiguity Resolution 
There are two types of ambiguities in the differential ranging system. 

15 First, there is ambiguity that is created if the signal is received by only three 
stations. In that case the hyperbolas are intjcrsecting in two points. 
Mathematically, that ambiguity cannot be resolved unless we have an a priori 
information (for instance - one of the points is in the sea), tracking data (the true 
intersection will move on roads) or other means such as Interferometer Direction 

20 Finding. The second ambiguity is the ambiguity in estimating AR. The reason for 
the ambiguity occurrence in AR is apparent from formula (9) which should be 
expressed as: 

AR,. = R, - R. = (A(p, - A(p. + 27in)* C/Aw (|0) 

25 

where n represents the fact that the difference in wavelength between the two 
frequencies could be larger than 1. Formulating it (10) at the first base station: 

Define - (0 - 2*71) ' 
30 R, = ( n„ + /2n) * X = ( n,, + M^, /2k) *(>, AX) l = R, / C 
A,,(t) = Ao * Cos {w( n,, 4- M^,. /2n) * /JC + <|)j 
w*A / C = (w - AwO*(?. + A?.) = 2n 
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A,, (I) = A,. * Cos I 2n n,, + M',, + ^\ 

A,.{t) = A„ * Cos ({w - Aw) t + <|)„) = A„ * Cos 1 2n n,, + M',, + ^„ J 
A(p.. = 27r (n,, - n,,) + M',, - M',, + (|) - (j>„ = Qtt (n,,-ii,0 + AM', + (|) - <}>„ 

5 For the second base station: 

R, = ( n,, + /2rr) * = ( n,, + /2ir) + A?v ): i = R, / C (12) 
A,,(t) = A„ * Cos { 27in,, + 4',, + ^\ 

A,,(t) = A„ * Cos ({w-Aw)t + = A„ * Cos {27in,, + 4^,, + ^„ J 
10 A(p, = 27t(n,,-n2,) + Tjj-M',., + (t)-(t)„=27t(n,,-n„) + AM^ + <j»-<|)„ 

By subtracting the two measured phase difference we deduce: 

A (0, - A <p, = 271 (n,, - n,, - n,, + n,, ) + A4', - AH\ = A R,, *Aw/C (13) 

15 

In order for AR,, to be non-ambiguous n„ - n,, - n,, + n,, has to be 
0. AR,,< D. In Figure 7, the length of one side of a triangle is longer than the 
: difference between the other two sides. Summarizing it: 

20 Assume Aq), - A<p, < n . . (14) 

27tAf * AR,, / C < 7t 
Af < C / (2* AR,2) 
Af < C /(2* D) 

25 For D practically limited by 50 Km, Af has to be smaller than 3 

KHz. We chose to use resolution of 2.5 KHz. The problem is that the accuracy of 
AR,, is inversely dependent upon Aw - the larger Aw is. the smaller the error in 
AR,,. 

30 Propagated Wave Form 

The propagated wave form has to resolve two contradicting issues. 
It has a resolution of 2.5 KHz and it has a maximal difference in frequency in 
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order u> reduce ranginii error. The proposed wave lorm sugyesis lhal ihc fn;quenc> 
dilTc ■ ■" • " ■ 

iSw: 



dilTerence Aw will be an inieuer mulliplicalion of ihe non amhitiuous frequency 



5 Aw. = i * ft w { 1 c i s N 1 . N = Number of hops. (15) 

The phase difference is incremenlinj. theoretical ly from liop to hop by fttp. not 
necessarily at the same order. 

10 5q) = 8\v*AR/C ^^^^^ 

Acp, = AWi * AR / C - i ♦ Sw * AR / C = i*5(p 

The evaluation of AR is done by computing ft<p which has no 
ambiguities. The proposed detection algorithm is as follows: 
15 1 . Create a complex vector of V, = {V„ Cos (Aq) . ). V„ Sin (A<p . )}. 

Those values are the quadrate output of the Aq) ^ calculation and Aq) , 
does not have to be calculated explicitly. 

2. Perform FFT ( ). 

3. Find the peak of the FFT and compute its exact "frequenc\ ." The 
2" "frequency" is a measure of the phase increments, and provides 

excellent estimation of 

The algorithm automatically iitjproves the signal to noise ratio by the 
factor of the number of hops. There is no dependency between the hop frequency 
and the frequency difference of the number of hops. The process is actually 

25 sampling in the frequency domain. The form of transformation from the frequency 
domain to the time domain is IFFT (Inverse Fast Fourier Transform) rather than 
FFT (Fast Fourier Transform). It has no effect on the results other than 
interpretation. The starting frequency does not have to be fiw or 0 but rather an>' 
frequency of our choice. It is equivalent to sampling in the time domain which 

30 does not start at T„ = 0 which has no effect on the result of the FFT other than a 
phase shift. 
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Adapialion of Siandard NexNct'^^ to DilTerenrial Ranuinu 
The overall system throughput could be reduced due to two impacts. 
First, the power of each bit is reduced by at least 2, due to the fact that the power 
is divided between 2 frequencies. In fact the power will be further reduced due to 

5 the non-linearily of the power amplifier which creates spurious frequencies that will 
consume pari of the power. Second, the transmission of two frequencies requires a 
new frequency plan. The frequency resolution is reduced to 2.5 KHz (from 7.5 
ICHz) and a special care must be taken to align the hopping with the transmission 
of the two frequencies. 

10 The creation of the two frequencies is done using a regular mixer 

which will multiply the local oscillator by a sine wave in Af/2. The result will be: 

Cos (27rf,.ol + Wo) * Cos ( 27t{ Af/2}t -f (|), + <|>^^ ) = ( 1 7) 

0.5*Cos (27r{f..o + Af/2}t + {<j),,, + 4,,,} + (j,^) 

15 + 0.5»Cqs(27i{f,o • Af/2}l + {ct),o ^ >ol - 4>..J 

^ ±7t/2 (BPSK) 

As is seen, the power is distributed between the tw^o signals. The 
phase modulation is applied through the differential modulating frequency. Both 
20 frequencies are phase modulated in inverse which has no significant since the 
demodulation is differential. The frequency difference is Af. 

Design of the Remote Mobile Unit 
In Figure 9. the modulating frequency is generated by the 

25 microprocessor 801 from a port called "TCMP" (time compare) which is capable of 
generating a square wave pulse train. The pulse train is modulated by 0' or 180^ 
using a digital XOR gate 802. The modulated square wave is passed through a 
low-pass (or band-pass) filter 803 which clips its harmonics and produces a sine 
wave which produces the two frequencies in the mixer. The ranging transmitter ol^ 

30 Figure 9 is designed to avoid carrier rejection and high order harmonics. The 
power amplifier performance is selected to be extremely linear. This results in 
avoidance of leakage of energy to other channels. 
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I he ll- circuit for ihc ranging system is shown in l-iuurc 10, Vhc 
local oscillator downconvcrls the IF signals from 70 ± 0.65 MHz lo 21.4 ± 0.01 S 
MM/, in 36 KH/ channels. Two separate conveners are used lo operaie in parallel. 
In order to ensure phase continuity between hops (ihe phase dilTerence beiween the 
5 two frequencies must remain between hops) we must base the double loop 
synthesizer on DDS. (Direct Digital Synthesizer). 

The process of finding the phase difference between the two 
frequencies is based on convolution of the two signals: 

10 Aq) = arg( S, ® S,*) (18) 

where arg(*) is defined as the argument of the complex variable (•) The process is 
identical to the process of interferometer DF. The convolution is done on a 
complete hop (45 msec) after the data extraction and the necessary phase reversal. 
15 The phase data of each hop is transferred to the central station. The process of 
differential phase comparison is done in the central station. The location 
processing includes differential ranging computation as well as a combination of 
differential ranging and direction finding. In operation, differential ranging 
requires that periodic calibration be performed in real time. The accuracy of the 
20 system will be in the order of 300 nanoseconds, which mean that liming error in 
the range of 50 nanoseconds must be compensated. The most cost effective way to 
do it is to set up several transmitters in predefined places and to calibrate the 
timing offset accordingly. Assuming a short term drift of 10 '*' (high performance 
oscillator), the time between calibrations will be: 
25 T = 50* 10'*' / 10-'" = 500 seconds (10) 

The present system supports remote activation of the calibration 
transmitters. 



30 Accuracy of DTP A 

The accuracy of the DTOA system is determined by comparing it to 
a regular triangulation direction measurement system. The geometric shape 
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rcpresenling the points thai have the properly of the same dilTerence time of arrival 
relative to a couple of bases is a hyperbola as shown in Figure 7. Accordintz to 
iriijionometric law of cosines: 

5 AR = R, -R, = {R- + (D/2))- + RD Cos (G)) ' - -{R- +{D/2)- - R D Cos(0) ] ' - 
(2R- + 2(D/2)- - AR- ) ' = 4 I ( R- + (D/2)- )- - { ( R D Cos (6) }- ] = 

8R- (0/2)= - 4R- AR- - 4(d/2)- AR- = 8R= (D/2)- - 4 1 ( R D Cos (G) |- 
Cos (G) = (AR/D) ( 1 + D- /4R-)'- for (AR/R) « I. (20) 

10 In order to simplify the evaluation of the inaccuracy of the location 

due to the inaccuracy of AR estimation, the inaccuracy of the asymptote is 
calculated. This is valid when D/2R«1 . Assuming D/2R«1, the asymptote 
f ormula is: 

15 Cos (ej = AR/D (21) 

The asymptote has an ambiguity of ±9 due to its being an even function. 
Estimating the error in G due to the error in AR is: 

20 <7(Ga) = o( AR) / D * SinOO (22) 

For an example, let E(AR) = 50m. D=20 Km, £(9^) is 0.14'. It 
could be shown that the hyperbola error is very close to the asymptote error. The 
error at the transmitter location could be evaluated as: 

25 

a(R) = Ra(G) = a( AR) / {(D/R) Sin(e)} (23) 

Note that the error in the estimation of the DTOA is multiplied by a 
geometric factor (1/ (D/R) Sin(9) ). The constellation of the two receiving bases 
30 and the transmitter could be scaled up or down and the effect on the errors is the 
same. The term (1/ (D/R) Sin (G) ) is recognized as GDOP. 
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(UX)I> = Sin(U) (24) 

where CiDOP = ucomelric delusion of propcriics. 

Hrror Sources in DilTcrential Raniiinu 
5 External Thermal Noise is one source of error in ihe preseni 

invention. Trom equation (9) follows equation (25): 

a(AR,: ) = a( Acp,. )* C/Aw (25) 

10 For signal to noise ratio (SNR) above 10 (IB the system uses the 

approximation (see Figure 11): 

a(A(p) = 1/ (2 SNR)' - : ..a ( Acp,. ) = 1/ (SNR)' " (26) 

15 The error in Aq> is dependent upon the algorithm. In order to simplify the 

process, in the alternative, the following is an algorithm for the extraction of the 
frequency out of the FFT. Find the peak of the FFT and its strongest neighbor. 
The value of these peaks are marked as P and P*. The following is the amplitude 
distribution is a sine wave in FFT with boxcar window (signal power equals A„ 

20 which is a function of the signal to noise ratio induced by Acp ). 

A - Ao * Sin (71 . 5q) . (27) 
N * Acp - 27ik + 5<; 



0 



Since Sin (5c) = Sin ( n - 8q) the system deduces: 

AVA = 8q/ ( 71 - 5q ) (28) 
5c = 7iA7(A + A^) 

In order to evaluate the accuracy of the estimation, a point is chosen 
which is the worst regarding the power distribution in the cells of the FFT, At this 
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point A = A\ Assume lhal the noise is while i.e. a A = a A' = \n. The error in ftc 



is: 



10 



a(5c) = 71 |aA{ A- + A^-)' -/(A + A')-' - !7i/2V2! * \aA/A] (29) 

A = Ao* Sin ( 7t/2 ) = A , *2/7i 

SNR.. ^ (A (/ ciA)" 

a( 5c )= {n/2}-/(2 * SNR, )' ^ 

Since the input noise is distributed over N cells of the FFT: 

SNR,= SNR.*N (30) 



SNRi is the input signal to noise ratio. Finally: 
a(AR,,) ^ a( A(p,, ) * C/Aw = ( 7t/8 ) * C / (5f*N^'"- {2 * SNR. (31) 
15 where C = 300.000,000 m/s 

. N = 53 

5f = 2500 Hz, 

SNRi = 10 dB =10 (integrating over a complete hop) 
20 And the final result of the estimation of AR,, is: 

a(AR,2) ^ 27 m (32) 

Multinath as a Source of Error 

25 The effect of multipath is very different from the effect of the noise 

in two ways. First, the multipath is not necessarily independent, thus the averaging 
that is done through the FFT will not provide an improvement of N^-. Second, the 
signal to reflection ratio is rarely in the order of 10 dB. Several simulations were 
done in order to evaluate the effect of the niultipath and the results showed that the 

30 frequency hopping scheme randomizes the effect of the multipath. 
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Dopplcr Shift as a Source oT Urror 
The Dopplcr clTccl is only on Ihc iVcqucncy dilTcrLMice 
dilTcrcncc \vc have: 

Ar,,„,, = v„,,,//. = v.„,,Aiyc' 

5 Where: 

Aluo,,, = Maximum Dopplcr shift 
V^,, = Maximum radial velocity (assume 40 m/s) 
Ar= Maximum frequency difference (up lo 1 MHz) 

10 Result in: 

Af.... = 0.13Hz 

Which is negligible since it causes a bias of 0.006 radians in each 
1 5 hop. which results in: 

E{AR^3J = E{ Aq) )* C/Aw = EJ(Ac/N)] * C/Aw = 13.6'm 
E(p) - Error in parameter estimation 

20 RMU Frequ ency Inac curacy as a Internar Source of Error 

The modulation of the dual frequency in the Remote Mobile Units 
(RMU) is done with a low cost oscillator which has a possible drift of 100 PPM. 
The error causes a direct error of 0.0001 in the differential ranging (sec 32) which 
is an error of several meters in the worst case. Timing is another source of internal 
25 error. The most difficult bias error to overcome is the timing error between 

receivers. It is obvious that even an atomic reference clock with an accuracy of 10" 

results in 50 nanosecond error within 1.5 hours. The system design is based on 
calibration transmitters in predetermined places. 

Internal Phase and Frequency Errors is another source of internal 
30 error. Differential frequency errors could be practically eliminated by using a 
TCXO (temperature compensated crystal oscillator) with I PPM accurac\'. Phase 
errors could not easily be eliminated internally. But. if the frequency generators 



In 1 Mil'/ 

(33) 
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.rc phase- co„,inuo„.,. ,hc phase error creates a differencial ra„.i,„ bias error ,sce 
cquauon T32, which has ,o ca.ihra.c conUnaousI, h, .he caUhra,- ,„ ,ra„s„,i„ers 
Nine the possible yroup delay at the II' crystal niters. 

Wave Form 

The transmitted wave form is a Hnai 

^ frequency transmission with a 
dillerencc of minimum of 200 KHz to 350 KU-7 «mM . • 

j<.Mz. with resolution of 2.5 KHz. 

The Two-Frequencv M kgntjalRano^^ 
'0 The following description is the complete mathematical analvsis of 

the two frequency differential ranging algorithm taking into account all known 
sources of error. Let / denote the true time for the /th hop. where / = 0 
corresponds to the beginning of the /th hop. ,n the following equations, a dot 
above the variable indicates its derivative. The transmitter frequency at the 
15 beginnmg of the ;th hop is given by 



20 



25 



f... =/,,, +&/, + 8 f, fi . \)T/ 

where 

/,,,, is the basic frequency of the /th hop; 
St/, is the transmitter frequency bias; 
6./, is the transmitter frequency drift; 
T is the lime interval between hops. 

The two instantaneous phases of the transmitter frequencies are 

= (f,.. + 0.55 f, t^ + <^,„ ^ „^ 

= a... -f„... >t ^ 0.56 f. t= <^,,„ . , „; 



(Tl) 



(T2^/; 
fT2/)y 



where 
30 • 



L,., is the frequency of the DSB modulation at the /th hop; 
is the initial carrier phase at the /th hop: 
„ is the initial phase of the DSB modulation at the /th hop; 
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t\ (I) is ihc phase noise of ihc iransmiucr oscillator. 
The insianiancous ranges al the nh hop are 

II) = i\n + \\„ [I + a - \yr\ + o.5«,„ |/ + fi . i^y-j: 

= \'\u + yj..,> (' - VT+ ().5«,„ a . ],^-T- I 
+ + a - 1)7]/ + o.5«,„ /- 

= '\. + V, , / + 0.5«, „ r. 

( I 

where 

k = 1 .2 is the receiver number; 

is the range between the transmitter and the kih receiver at the 
beginning of the first hop; 

v,„ is the velocity along LOS between the transmitter and the Anh receiver 
at the beginning of the first hop: 

«,„ is the acceleration along LOS between the transmitter and the fah 
15 receiver at the beginning of the first hop: 

/•,, is the range between the transmitter and the Ath receiver at the 

beginning of the /th hop; 

• 'V, is the velocity along LOS between the transmitter and the Anh receiver 

at the beginning of the /th hop. 

20 The two received signals at the Mh receiver are 

,.. <V = 2-' - exp(/27t4,.„ 0 - K^;(l)))m,„ (/; + (,) fT4a) 
(I) = 2-' - exp{,-2K(l>,, (t - K,r,(imm„. (i) ^ u-,, (n . (T4h) 

25 where 

/c„ is the inverse of the speed of light: 
m,„W and n^O are the muUipaths of the two iVequcncies: 
u,,, fO and m;,(0 are additive white noises for the two frequencies. 



30 



In practice, the two frequencies are not received separatelv. but their 
sum „ 0) -H (0 is received. However, since the two frequencies are separated 
b> 7.5 KHz at least, which is much above the bandwidth of each frequency, we 
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assume that they can be separated by prelllieriny. This also explains ihc Tacior 2' 

namely, lo rcllcel the facl that ihe lolal iransmiuer power is divided beiween ihe 
iwi) frequencies. 

The kih receiver generates two reference signals whose instantaneous 

5 phases arc 

<t>r,. (0 = (f\o^L,JO -6 ^ b t,0 4),,,, (T5^/) 
<i>rkH (0 = (fiu'LjU -0 + 6 t,,) + 4>,,„ , (T5b) 

10 where 

6 is the delay, with respect to the beginning of the /th hop, at which 
demodulation starts; 

6 t, 1^ is the time bias of the k\h receiver. 

<j),j^ (, is the initial phase of the reference signal of /ih receiver at the nh hop. 
15 Note that the receivers are assumed lo have "ideal" frequencies (zero 

frequency deviations), because the reference signals are synthesized digitally, so 
and /;„ , are just numbers in the computer. 

The demodulated signals at the receivers are 

20 cl,^, (V = exp {/27i(|>_ 0) 1 . (t) (T6a) 

d,, (!) = exp {/2n(t),,, (0) \, (i) . {TGh) 

The equations derived so far are used in the simulation. In order to 
derive the estimation let us now ignore the noise. Also express the multipath in 
25 terms of its magnitude and phase, that is. 

0) = A,^, (I) exp {/27iv|i,„(t); (T7^;) 

niu. (t) = A,, fi) exp {/2n\i/,,(t);. (T7h) 

Then, 

30 d,, (0 = 2 '- A,,, (i) exp {jln (i) . 4,^^^ a . ^.r, dj) ^ y^^JO]) (T8a) 

ciu. (0 = 2-' - A,, (i) exp [fin [^^,, (f) . <i>_^, (i - k,,/*, (dj - 4',^^)]]^ (T8/» 
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The insianlaneous phase oT each demodulaicd sitznal contains a 
quadratic icrni. a linear term, a constant term, and a random term. Denote the 
cociricicnts of the linear terms by ./,,,, and/,,, respectively, h is easy to verity thai 
the linear terms are approximately given by 

5 

.U. = K„ v„ (U, + ) - 6/;- 6 i; (i - 1 )7- (T9«) 
U = xo (L. - .L.. ) - ¥r ft i; (i - 1 )T. (T9b) 

Suppose we estimate the linear terms (using a Discrete Fourier 
10 Transform with interpolation), and use the estimates to form 

/s A /V A 

a, =./di„ - L2.. * K«X.,rv,, - V,, ) « K„./„/v, , - V,, ). (TIO) 

The use of a , is explained below. 
15 Next let us form the signal 

c/> = </*.,r'M** = (t) A,at) 

exp{/27r [ct>,H 0 - k„ r, (t)) - <t>,„ (I • k„ r, (t)) + <i),,„ (t) - <t>ivy, (I) + 

20 

We have 

<^),, r/ - (0) - <i>.. 0 - K„ f, (0) = -2f„,. 0 ■ K„ r, 0/} - 2<t>,.,, 

= 2/,,,, K„ t\. - 2j,„, (1 -K„ V,, )t + /;„>•„ r - 2<t>„,, (Tl 2) 

25 <I>H.. (0 - <t>ri, (0 = 2L, t - 2f„,,, e + 2f„,, 6/,, . (Tl 3) 



So. 



(0 = 0.5/lt„ (t)A,, (I) 
30 exp{/2Tr {2f,,, .K,r^ ^%L ,>ioVkjH.y,'<-,P 5 • 



(T14) 
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The maximum value oi" 2f„,,K„\\ / is aboul ?-\(y\ and ihc maximum value or/,„,K„ 
a^j- is aboul 2-10 \ so ihcsc two icrms can be neulceled. ieavinii 

5 

Assuming that /n^,, (f) and /'z; are strongly correlated, \vc have i/;., ^ xj^u^. Let 
[/^^, 0 < /7 < A' - 1 } be the sampling instants of the /th hop and form 

(T16) 



p. y= (2n) -^-angle 



ln=0 



10 and 

P. = iS,, - * 2/„,,K„Vr,, - /•,;,) + 2/:,,,(6/h - 6/,:)- (T17) 

The bandwidth of the signals d^^/tj and €1^^0)13 limited to about 200 
15 Hz (twice the maximum Doppler frequency). Therefore they should be low-pass 
filtered to this value and decimated to a sampling frequency of aboul 20D'^Ife 
With the present value of T (50 milliseconds), we gel 10 samples per hop. If the 
bit rale is also 10 per second, the range estimation algorithm can be conveniently 
combined with the bit-detection algorithm, since the phase of each bit is needed for 
20 both purposes. The simulation does not include the filtering, but the standard 
deviation of the additive noise is computed in agreement with the assumed 
bandwidth, 

Finallv use the estimates a, from equation TIG to compute 
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/-I 

V = w - y:.fL "1 

1 

10 As we see, the role of a, is to compensate for the motion during the 

hops and to estimate the initial range difference. Motion compensation depends on 
the receivers carrier frequencies being identical. If they cannot be made identical, 
this compensation will not be performed, and will be used for subsequent 
processing instead of y,. 

15 The final stage of the algorithm is the estimation of the range 

difference from {y,. 0 < / < I - 1 }, (or from {P,, 0 < / < 1 - 1 {, as discussed 
below). Let be the ordering of the hops such that (i) =X„„ + .0 < f 

< / -1. Let Ti f = expt/27iy . Then, except for the receiver synchronization 
errors b t . r\ ^ \s given by 

20 

Ti , = exp(/47i (r,o - /So) LLn + (T19) 

Therefore. 2/^,,j/c,Xr,„ - /s.o) can be estimated by performing a Discrete Fourier 
Transform with interpolation on the sequence {r| ^, 0 < f < / - 1 }. 

25 If we use jS, instead of 7,, the sequence of ri, will not be a pure 

sinusoid (even in the ideal case), due to the transmitter motion. However, if f(i) 
is sufficiently random, the phase deviation seen by (using a Discrete Fourier 
Transform with interpolation) will be random, so this should not have much effect 
on the estimate. For example, if the speed is 30 m/s, the transmitter motion during 

30 the entire sequence of hops is no more than 120 m. With /„,„, (/ - 1) = 60 

KHz, this translates to about 0.05 •27c radians phase error. It is therefore expected 
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that ihc Discreic I-ouricr rranslbrm with inicrpolation procedure will eivc a good 
cslimale of ihc avera^^c range, or ihc range al Ihe (//2)th hop. 

The MuUinath Model 
5 Mullipalh is basically a linear phenomenon, lhal is. the received 

signal is a supcrposilion of several signals, corresponding lo several propagation 
paths from the transmitter to the receiver. One of these is the direct path, whose 
impulse response is assumed to be an impulse, with a gain factor a and a zero 
delay (the nominal propagation delay lo the receiver is immaterial, and can be 
10 ignored). All the rest are reflected paths. Each reflected path is characterized by a 
random gain and a random delay. These random parameters vary over time in 
general, so the total response at time / to an impulse transmitted at lime /(, is given 
by 

h (t. t,) =a6 (t't,) -^ccn(^o) <5 (t:-t,'T,(to) ) . ( T20) 



15 The case = 0 is called Rayleigh faclinfi^ w'hile the case ^/ > 0 is called Rice fadin}-. 
In the Rayleigh case there is no direct propagation, only reflections. This happens 
when there is no line of sight between the transmitter and the receiver, and the 
wavelength is loo short for diffraction to have an effect. 

When the transmitted signal is s(t). the received signal is given by 

20 the convolution of .s//; with h(l. that is. 

z(t) =as(t)^Ea,(t)s(t'V,(t)} . iT21) 
n 



If we now assume the existence of infinitely many reflected paths, each having an 
infinitesimal contribution. (T21) becomes 



z (t) =as (t) Hoc (x, t) s (t'X)dx, 

J 0 



{T22) 



25 where a(Tj) is the infinitesimal contribution at time / of the refiection delayed by 
T. The statistical properties of afrJ) will be discussed later. 
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Now Icl be a C\V signaK namely. s(f) ="cxp( ;/27T/n^ i . Then 

J- (V) -as (t) *f^^ C/exp ( t -r) } gt: ^cxp * i 27i r\.t } {a-mit) ] , ( T23) 

where 

m(t)= raix, t;expl-j27ifo^}dt . {T24\ 

Jo 

As we see from (T24). Ihe effect of muUipaih on a C\V signal is 
nwldplicative, and il depends on the function of a(7j). Here we are interested in 
5 short-term multipath in urban areas and mobile transmitters (or receivers). A 

common model for afrj) in such scenarios is as a Gaussian process, stationary and 
narrowband in /. and nonstationary but uncorrelated in r. Thus we assume 
Ea(r,f) = 0, 

r(T, , X.: /, . /,) = Ecc(r, ,/,)a(T, . /,) = V, (/, . /.) (t,)5(t, - x,). (T25) 
10 The temporal correlation r, - depends mainly on the physical 

motion of the transmitter. In particular, if the doppler frequency of the transmitter 
is f,i„f,fyn.r (corresponding to the wavelength and the transmitter speed), then the 
spectrar density of T, 0^ - is limited to ±/,„,,,„,, . The shape of the spectral 
density is difficult to estimate, since it depiends on the given scenarW: mainly on 
15 the spatial distribution of the reflections. In the simulation we used a flat spectrum 
in the range ±Lf,pf,r simplicity and "universality." 

The function (xi) describes the relative multipath power at 
different delays. It is also scenario-dependent and difficult to estimate. However, 
based on the graphs known in the art, we assume a functional dependence of the 
20 form 

(x,) = exp{.x, /x^/,, (T26) 
where t„,^ is a scenario-dependent parameter. Note lhal when T. rr,) is plotlcd in 
dB, il shows a decreasing straight line. 

It is interesting to compute the multipath correlation for both 
25 different times and different frequencies. Let w(/, . /,) be the multipath at 
frequency /, and time . and similarly for , /.). Then 



wo 96/25673 PCT/GB96/00270 

41 



r„, ( t, , t_..; ,f^)=Em{t^. f. ) m* { t.. , f .) 



10 ° 



1+J2TI (/j-fa) T 



inp 
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10 



Table 1 shows the results of simulation of the present system with 
the given parameters. The frequency is selected as described above as/„„. 
(/ - 1), SNR is the signal to noise ratio. MP = 0 indicates Rayleigh fading 
15 and MP = 1 indicates Rice Fading. STD is the standard deviation (one a), and the 
resolving power of the algorithm is given as minimums (Min) and maximums 
(Max). 



20 Although specific embodiments have been illustrated and described herein, 

ii will be appreciated by those of ordinary skill in the art that any arrangement 
which is calculated to achieve the same purpose may be substituted for the specific 
embodiment shown. This application is intended to cover any adaptations or 
variations of the present invention. 
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CLAIMS 

1 A mobile ininsinilier localion syslcm. coniprisiivj: 

i\ luiibilo iratisiuillor havinu a uniqiio idcnlillcr. ami lia\ inii ;i 

iransmillcr lor sitnuliancDusly soiulinti a first frcqucncy-hoppcd spread speclriim 
5 carrier signal and a second rrcqucncy-hopped spread speelrum carrier signal, ai 

least one of which containing tlie unique identifier, the first frequency-hopped 

spread spectrum carrier signal having a different frequency than the second 

frequency-hopped spread spectrum carrier signal; 

a first base station having an antenna operable for receiving the first 
10 and second frequency-hopped spread spectrum carrier signals and capable of 

determining a first lime of arrival of the frequency-hopped spread spectrum carrier 

signals; 

a second base station also having an antenna operable for receiving 
the first and second frequency-hopped spread spectrum carrier signals and capable 
15 of determining a second time of arrival of the frequency-hopped spread spectrum 
carrier signals; and 

a central station connected to the first base station and the second 
base station, and operable for: 

1. ) receiving the first time of arrival of the first and second 
20 frequency-hopped spread spectrum carrier signals from the first base 

station: 

2. ) receiving the second time of arrival of the first and 
second frequency-hopped spread spectrum carrier signals from the 
second base station; 

25 3.) generating a first hyperbolic plot based on the differential 

time of arrival of the first and second frequency-hopped spread 
spectrum carrier signals at the first and second base stations; and 

4.) using a third information source to determine the location 
of the transmitter along the first hyperbolic plot. 

30 

2. The mobile location system according to claim 1 wherein a reference station 
operates periodically to transmit a reference first and second frequency-hopped 
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spread spectrum carrier signal from a known location to calibrate the first and 
second base stations and the central station. 

3. The mobile location system according to claim 1 wherein the first base 
5 station and the second base station both include atomic clocks to synchronize the 
time of arrival calculations. 



4. The mobile location system according to claim 1 wherein the third 
information source is a geographic map of the area around the base stations which 

10 provides information on the probable location of the transmitter along the first 
hyperbolic plot. 

5. The mobile location system according to claim 1 wherein the third 
information source is a third base station which calculates the time of arrival of the 

1 5 first and second frequency-hopped spread spectrum carrier signals and wherein the 
base station generates a second hyperbolic plot to intersect the first hyperbolic plot 
to identify the location of the transmitter. 

6. The mobile location system according to claim 1 wherein the transmitter 
20 transmits a plurality of hops each having a different first and second frequency- 
hopped spread spectrum carrier signals and the time of arrival is determined by 

a.) calculating the phase difference between the first and second frequency- 
hopped spread spectrum carrier signals for each of the plurality of hops to produce 
a plurality of phase differentials: and 
25 b.) performing an Inverse Fourier Transform on the plurality of phase 

differentials to produces a time of arrival. 

7. A method of locating a mobile transmitter location, comprising the 

steps of: 

30 receiving at a first base station a dual radio carrier having two 

frequencies and calculating therefrom a first time of arrival of the dual carrier; 
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receiving at a second base station a dual radio 
carrier having two frequencies and calculating therefrom 
a second time of arrival of the dual carrier; 

comparing the first time of arrival with the second 
5 time of arrival; and 

determining the location of the origin of the dual 
radio carrier by hyperbolic intersection location • 

8. A mobile transmitter having a unique identifier for 
10 use in the method of claim 7 comprising a transmitter for 
simultaneously sending a first frequency-hopped spread 
spectrum carrier signal and a second frequency-hopped 
spread spectrum carrier signal, at least one of which 
containing the unique identifier, the first frequency- 
15 hopped spread spectrum carrier signal having a different 
frequency than the second frequency-hopped spread 
spectrum carrier signal . 

9 . A base station for use in the method of claim 7 
20 comprising an antenna for receiving the dual radio 

carrier and means for determining the first or second 
time of arrival of the dual carrier. 

10. A mobile transmitter location system comprising: 
25 a first base station having an antenna for receiving 

a dual radio carrier having two frequencies and means to 
calculate therefrom a first time of arrival of the dual 
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carrier ; 

a second base station having an antenna for 
receiving a dual radio carrier having two frequencies and 
means to calculate therefrom a second time of arrival of 
5 the dual carrier; 

comparing means for comparing the first time of 
arrival with the second time of arrival; and 

location determination means for determining the 
location of the dual radio carrier by hyperbolic 
10 intersection location. 
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